Seeing topological order in time-of-flight measurements 



E. Alba,^ X. Fernandez- Gonzalvo,^ J. Mur-Petit/ J. K. Pachos,^ and J. J. Garcia-Ripoll^ 

^Instituto de Fisica Fundamental, IFF-CSIC, Calle Serrano 113b, Madrid 28006, Spain 
^School of Physics and Astronomy, University of Leeds, Leeds, LS2 9JT, United Kingdom 

(Dated: January 20, 2013) 

In this work we provide a general methodology to directly measure topological order in cold 
atom systems. As an application we propose the realisation of a characteristic topological model, 
introduced by Haldane, using optical lattices loaded with fermionic atoms in two internal states. We 
demonstrate that time-of-flight measurements directly reveal the topological order of the system in 
the form of momentum space skyrmions. 
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Different phases of matter can be distinguished by their 
symmetries. This information is usually captured by lo- 
cally measurable order parameters that summarize the 
essential properties of the phase. Topological insulators 
are materials with symmetries that depend on the topol- 
ogy of the energy eigenstates of the system [1 . These 
materials are of interest because they give rise to robust 
spin transport effects with potential applications ranging 
from sensitive detectors to quantum computation [2]|3j. 
However, direct observation and measurement of topo- 
logical order has been up to now impossible due to its 
non-local character. Instead, experiments have relied so 
far on indirect manifestations of this order, such as edge 
states and the quantization of conductivity. 

Ultracold atoms facilitate the implementation of arti- 
ficial gauge fields [4 . Here we distinguish proposals that 
generate continuous fields [5 , such as the recent exper- 
iment by Lin et al [6 , from those that rely on optical 
lattices and engineering of hopping [7 . We will concen- 
trate on the latter, introducing a method based on stan- 
dard time-of-flight (TOF) measurements that can iden- 
tify a topological character in the quantum state of the 
system. Our starting point is a possible implementation 
of Haldane's model using fermionic atoms in two inter- 
nal states. The topological nature of its ground state is 
witnessed by the Chern number. This number counts the 
times the ground state, written as a spinor, wraps around 
the sphere, as a function of momentum. We demonstrate 
that TOF measurements reconstruct the Chern number 
in a way which is robust against the presence of external 
perturbations or state preparation. Our method can be 
adapted to other quantum simulations of topological or- 
der in optical lattices jSHTT] , as many already use internal 
degrees of freedom of the atoms to encode the order. 

One common mechanism for the appearance of topo- 
logical order is based on the topology of the eigenstate 
manifolds. Consider a real-space lattice whose unit cell 
has d quantum degrees of freedom — position of the par- 
ticle, spin, etc — . Its energy band description, £^^(k), 
has d eigenstates, '?Am=i...(i(k), per value of momentum 
k = (kx-,ky) in the Brillouin zone, B. Different config- 
urations of the vector fields = z(?/^^(k)| V/eV^rn(k)), 



defined as the gradient in momentum space of the wave- 
functions iprni are characterized by topological invariants 
such as the Chern number of each band ^, 

1^m = ^ I ^kX ^md?k. (1) 

In the case of the quantum Hall effect, the energy bands 
are separated from each other and the material becomes 
an insulator for appropriate Fermi energies, Ep. In a 
real setup, with finite boundaries, the sample can have a 
quantized non-zero conductivity given by the topological 
invariant = I^^e^<Ef which is a signature 
itself of topological order [18 . The transport is then sup- 
ported by "edge" states that are localized on the bound- 
ary of the material, with interesting properties, such as 
spin polarization or chirality, insulation from noise and 
resilience to perturbations [3 . 

There are two essential routes towards topological or- 
der in momentum space, depending on how we realize the 
d quantum degrees of freedom mentioned above. One 
is to start from charge carriers with intrinsic angular 
momentum and a spin-orbit coupling, as in the Kane 
model [1 or in semiconductor structures [19]. A differ- 
ent approach is exemplified by the Haldane model [20] , 
built on a honeycomb lattice where the unit cell has two 
sites and a spatially modulated magnetic field breaks the 
parity. We shall propose a generalization of Haldane's 
model that relies on optical lattice technology and two 
atomic hyperfine levels (pseudospin) to distinguish be- 
tween the sites of the unit cell. This enables us to directly 
extract the Chern number from the spin textures in TOF 
images [21] and demonstrate the topological order. 

Let us consider a honeycomb lattice constructed out of 
two triangular sublattices, A and B [Fig.jl^]. Each of the 
sublattices hosts fermionic alkali atoms in a different in- 
ternal state, I a) and The model is parameterized by 
four couplings: the hopping amplitudes inside the same 
species lattice, ta and t^, the energy difference between 
A and B sublattices, £, and the coupling between sublat- 
tices, tj/c, which can be induced by a Raman laser and 
controlled at will [22]. If the lattice is deep enough and 
the tunneling amplitudes and interaction energies remain 
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FIG. 1: Haldane-type model, (a) Two triangular optical lattices (A and B) are Raman coupled by a laser that makes an atom 
switch sublattice {A ^ B). This allows for both next-nearest neighbor hoppings, ta,b, and a complex nearest neighbor hopping, 
tjkj whose phase depends on the momentum imparted by the Raman laser, Ap. (b) Phase diagram of the zero-energy ground 
states, as a function of the energy imbalance between lattices, e/ta, and the momentum imparted by the laser, Ap = (0, Apy). 
We plot the exact phase boundary in the thermodynamic limit (black solid line), together with a color graded simulation of the 
Chern number for a finite lattice with 2 x 17 x 17 sites and 20 x 20 pixels. The diagrams on the right hand side show how the 
Dirac points are displaced on the distribution of Bloch vectors Sz(k) induced by the Hamiltonian (blue negative, red positive). 



small compared to the interband separation, we may use 
single band tight-binding model [33 

H= {Ujblaj + t^aj6,) + ^ s{a]aj - b]bj) 

{ai,bj) j 

In the presence only of the hopping tij the energy spec- 
trum consists of two energy bands that meet at two 
"Dirac points". At half filling the low energy physics of 
the system is dominated by the linear dispersion around 
these points, the "Dirac cones" [23]. Due to the pres- 
ence of ta^h and e, the effective Dirac fermions acquire 
a mass that depends weakly on momentum. The posi- 
tion of the Dirac points on the mass landscape determines 
whether the model is topologically ordered or not. In our 
cold atoms simulation this is controlled using the Raman 
lasers to attach a phase to the hopping [7 

tjk ^ tex.p{i(l)jk), (l)jk = Ap • (x_^- +x/e)/2. (3) 

This phase displaces the energy bands created by the Ra- 
man hopping, relative to mass landscape generated by 
the other contributions, ta^b and £, as shown in Fig. 
When the Dirac points have opposite signs of the mass, 
the Chern number ([T]) automatically becomes non-zero. 
Intuitively, while the total flux over each hexagonal pla- 
quette is zero, the bipartite nature of the lattice allows 
the phases (j)jk to have a non-trivial effect: along the path 
1 ^ 2 ^ 3 ^ 1, depicted in Fig. [l^ the local effective 
magnetic flux, 0i2 + 023, is also different from zero. 

The momentum space Hamiltonian associated to ([2| 
has the structure given in the introduction ^ 

H(k) = x(k)l-^(k)S(k)-^T, (4) 



with two energy bands, ±£^(k), the Pauli ma- 
trices a = (a^^a^^a^) and a normalized pseu- 
dospin S(k) labeling the state of the atoms 
in the {|a), \b)} space. In our model S(k) oc 
(tRe/(k-Ap),tIm/(k-Ap),£+(ta-t5)^(k)), with the 
complex functions /(k) = Xln^o i 2 e~*^"^'^^ and ^(k) = 
^^^3 4 5 cos(k • v^a), a set of displacements G ^ x 

{(-2, 0), (1, V3), (1, - V3), (0, 2^3), (3, ^3), (3, - ^3)} 
and the honeycomb lattice spacing a. The energy shift 
X = (^o + h)9 does not affect the topological phase. 

The topological properties of the model can be ob- 
tained from the field S(k). In particular, the lowest en- 
ergy band has a total Chern number 

iy=^ [ S-{dk^Sxdk^S)d^k. (5) 

Fig. summarizes the three different phases that can 
be accessed by means of the effective magnetic flux, 
= (O5 ^Py)^ and the imbalance between lattices, e/ta. 
First of all we find a trivial region, u = 0^ which is topo- 
logically equivalent to graphene with a mass term. When 
we interpret the associated spin texture as a map onto the 
Bloch sphere, both cones have the same effective Dirac 
mass and point to the same pole, Sz > 0. Since they cover 
the same polar cap in opposite senses, u = ±{^ — ^) = 0. 
Moving across the solid black line in Fig. [l]3, the lattice 
undergoes a quantum phase transition into a topolog- 
ically non- trivial phase. Now the cones on inequivalent 
Dirac points are positioned at opposite poles of the Bloch 
sphere, forming a skyrmion [18 [Fig. that covers the 
whole sphere and thus = + ^) = ±1. 

The setup in Fig. [l^ may be experimentally realized 
along the lines of Ref. [24 , combining spin-dependent 
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FIG. 2: (Color online) (a) Spin texture of the Hal- 
dane model, interpreted as a mapping from momentum 
space, {kx,ky), onto the Bloch sphere, S = (cr) oc 
(cos sin ^, sin sin ^, cos The colors and arrows show the 
polar angle ^(k), and the azimutal component of the spin, 
(Sx, Sy), for a topological phase u — (b) Interference pic- 
tures for this phase which result in momentum density distri- 
butions for the a (up) and b (down) particles. A square lattice 
of 20 X 20 pixels partitioned into triangles (yellow) is used to 
compute the estimate vd — 0.9 [Eq. The enclosed area 

corresponds to (a). 



potentials [25l [26], with recent techniques for creating 
dipole traps using microscope objectives [27]. We sug- 
gest projecting two triangular lattice patterns on a two- 
dimensional sheet of light that traps the fermionic atoms. 
An electro-optic phase modulator controls the relative 
displacement of the lattices [28j and the appropriate 
weights of left and right circularly polarized light [25l [26] . 
The result is two hyperfine ground states of the same 
fermionic species confined on the two triangular sublat- 
tices of the honeycomb pattern. Thanks to this con- 
figuration, the distribution S(k) oc (cr) can be experi- 
mentally determined from the TOF images that appear 
when the atoms are released from the optical trap. A 
typical experiment would begin with a Mott state in 
which only the A sublattice is filled, and adiabatically 
progress to larger values of ta-, t and e. Once the ap- 
proximate ground state is prepared, switching off the 
trap in adequate timescales ^ projects the atom cloud 
into the momentum density distributions, 71^,5 (k), giv- 
ing direct access to one of the pseudospin components 
^'^(k) = ^[na(k) — 77,5 (k)]/[na(k) +715 (k)]. A fast Raman 
pulse during TOF allows us to rotate the atomic states 
and map Sx and Sy to 6'^, reconstructing the whole vec- 
tor field. Actual experiments "pixelize" the time of flight 
images, counting the number of atoms on each "square" 
of the effective Brillouin zone and estimating the averages 
of Sx^ Sy or Sz' Either through repetitions or through 
self- averaging in an experiment with multiple copies of 



the lattice, we will obtain a set of normalized vectors 
{'^m}^=i, evenly sampled over momentum space. As 
shown in Fig. |2]3, we suggest identifying the pixels with 
the nodes of a triangular lattice, T = {S^^, S/e^, S/^}, 
approximating the integral v by its discretization 

i^D~^Y. ■ Sfe. y^Si^ = v + 0(47rVL2). (6) 

T 

The value vd has the properties of a topological quantity 
— stability and robustness against local perturbations — , 
and is also stable with respect to the discretization [29] , 

We have compared the thermodynamic limit distribu- 
tion S(k) with realistic finite-size lattices with imperfec- 
tions. For this we have exactly diagonalized Eq. ([2| on 
a finite lattice with up to 20000 sites, including the ad- 
ditional harmonic confinement term, ^Tncj^xf, which is 
typical from cold atom experiments. Our plots report 
simulations with t/h = 1 kHz, ta/t = 0.5 and tb = 0, 
using r = muj'^o? /2t to parameterize the influence of 
the harmonic confinement. Realistic values for a lattice 
with a ~ 400nm range from r = 10~^ (^Li in a trap 
with uj/2iT = 60 Hz) to r = 0.02 ("^^K in a trap with 
UJ/27T = 100 Hz), but we probed up to r = 0.25. The 
results are very insensitive to the number of atoms, as 
already for 17 x 17 sites the interference pattern provides 
the right phase diagram [Fig. ^p]. The Chern number 
is also very robust with respect to the discretization: a 
20 X 20 pixelization deviates from the theoretical value of 
u only 10% [Fig. [3]d], in line with the error 0(47rV^^), 
expected from a discretization with a smooth integrand. 
Moreover, still captures the discontinuity across the 
topological phase transition [Fig. ^p]- Contrary to the 
global density of states and the eigenenergies, the ap- 
proximate Chern number ujj is also robust against inho- 
mogeneities. Already for a confining trap with r = 0.02 
the Dirac cones are no longer evident [Fig. [3^] , but the 
Chern number is still close to ±1 with a good signal-to- 
noise ratio [Fig. [sj^]. This is due to the wedding cake 
structure introduced by the harmonic trap [30] : for suffi- 
ciently strong traps there is always one ring or disc host- 
ing ntop atoms in a topological phase [Fig. [sj^-d] . Only 
these atoms contribute to the total Chern number, much 
like only superfluid atoms add to the interference peaks 
in experiments with bosons in optical lattices [ST. Note 
also how, as shown in Fig. [Sji, for low densities there are 
not enough atoms to form a topological phase and the 
Chern number deviates from ±1. However, raising the 
trap brings the chemical potential up to a level in which 
the ffist disc with particles in a topological phase is cre- 
ated, and converges to ±1. Finally, we expect also 
a good behaviour in finite-temperature simulations be- 
cause the effect of temperature only changes the length 
of vector {a) (i.e. the signal), but not its orientation S. 

Summing up, we have presented a robust and very gen- 
eral method to detect topological order in momentum 
space using ultracold atoms in various internal states and 
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FIG. 3: (Color online), (a) Density of states and (b) Chern number simulation at ejia — —0.5 for a lattice with r = (black), 
10~^ (dashed) and 0.02 (dotted), (c) Chern number for ejia — —0.5 and 'Py — Z'n j^a (blue, above) or 27r/a (black, below) vs. 
confining potential strength, r, starting with 1/2 (solid) or 1/4 particles per site (dashed) on a lattice with 50 x 50 sites. The 
values obtained using 50 x 50 pixels are compared with those from a 20 x 20 matrix (gray), (d) Density plots for the points 
marked in (c), showing the wedding cake structure where regions with n = 1 do not contribute to the Chern number. 



TOF images. As a very relevant application we have in- 
troduced an experimental proposal to generalize the Hal- 
dane model [20|. We found that the topological phases 
and our method are both robust under (i) use of small 
finite lattices (ii) coarse grain measurements of the spin 
texture, (iii) inhomogeneous potentials superimposed on 
top of the lattice, and (iv) errors in the exact values of 
the chemical potential, number of atoms or finite temper- 
ature. We believe our proposal is thus advantageous with 
respect to other indirect detection schemes — edge trans- 
port, eigenstate preparation [TT] or local estimates of the 
density of states ^30^ — , which may be more sensitive to 
temperature and imperfections. Compared also to the 
numerical protocol in Ref . [32] , our method only requires 
a single set of measurements instead of reconstructing 
spatially dependent correlators in position space, which 
are not easily accessible in optical lattices. Finally, the 
implementation of our ideas would represent the first di- 
rect visualization of non-local topological order. 
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